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Abstract We report on variations in the mean posi- 
tion angle of the 20 milhsecond pulsars being observed 
as part of the Parkes Pulsar Timing Array (PPTA) 
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project. It is found that the observed variations are 
dominated by changes in the Faraday rotation occur- 
ring in the Earth's ionosphere. Two ionospheric models 
are used to correct for the ionospheric contribution and 
it is found that one based on the International Ref- 
erence Ionosphere gave the best results. Little or no 
significant long-term variation in interstellar RM was 
found with limits typically about 0.1 rad m~^ yr~^ in 
absolute value. In a few cases, apparently significant 
RM variations over timescales of a few 100 days or more 
were seen. These are unlikely to be due to localised 
magnetised regions crossing the line of sight since the 
implied magnetic fields are too high. Most probably 
they are statistical fluctuations due to random spatial 
and temporal variations in the interstellar electron den- 
sity and magnetic field along the line of sight. 
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1 Introduction 



ISM: general — radio 



The observed radiation from pulsars is affected by prop- 
agation effects along the path between the pulsar and 
the observatory. One of the various propagation effects 
is Faraday rotation, that is, the rotation of the plane of 
linear polarization caused by the different phase veloc- 
ities of the two hands of circular polarization in a mag- 
netised plasma. Pulsar radiation typically has strong 
linear polarization making it relatively easy to observe 
Faraday rotation for these sources. Faraday rotation is 
described by the rotation measure (RM), defined by 



= RM 



(1) 



where ^ is the linear polarization position angle (PA) 
and A = cjv is the radio wavelength corresponding to 



2 



radio frequency v. The rotation measure is approxi- 
mated by 



RM = 0.810 / ne-B ■ dl 



(2) 



where is the interstellar electron density in units of 
cm~^, B is the vector magnetic field in microgauss, dl 
is an elemental vector along the line-of-sight toward us 
(positive RMs correspond to fields directed towards us) 
in parsecs and RM is in units of rad m~^. 

Along with the dispersion measure (DM) , defined by 



DM 



D 



Up dL 



(3) 



where DM is in units of cm~'^ pc, the RM may 
be used to determine the strength and direction of 
the i r iterstellar magn e tic field (see, e.g., Manchester 
19741 : iHan eraD bood iNoutsos et al.]l2008l ). Combin- 
ing Equations [2] and El we obtain an estimate of the 
mean line-of-sight component of the interstellar mag- 
netic field weighted by the local electron density along 
the path: 



= 1.232 — uG, 

\ 11/ Y)M ^ ' 



(4) 



where RM and DM are in their usual units. 

Long-term changes in pulsar PAs can result from 
changes in the polarization of the emitted radiation or 
from changes in the RM along the path (Equation ([T])). 
Intrinsic PA changes could result from precession of the 
pulsar spin axis or changes in the magnetic configura- 
tion of the star. RM changes can occur as the path 
to the pulsar traverses different regions of the interstel- 
lar medium (ISM) or in the Earth's ionosphere due to 
the diurnal and other changes in the ionospheric total 
electron content. 

Correlated RM and DM variations with time have 
been observ ed for several pulsars. 



Hamilton et al 



(|1977L Il985r ) analyze the variations of RM and DM for 
the Vela pulsar, finding an approximately linear vari- 
ation in both with slopes of 0.73 rad m^^ yr^^ and 
—0.04 cm~^pc yr~^ for RM and DM respectively. They 
interpret the RM and DM variations as the result of the 
relative motion of a magnetized filament of the Vela su- 
pernova remnant across the line-of-sight and estimated 
the electron density and the effective mean line-of-sight 
component of magnetic field in the filament. Adopt- 
ing 500 km s~^ as the transverse velocity of the source, 
they obtain a transverse scale length of 5.1 x lO"** pc 
per year. Assuming that the line-of-sight dimension is 
comparable with the transverse scale gives an electron 
density in the region of Arig = 78 cm~^. The mean 



line-of-sight component of magnetic field in the varying 
region is given by 



= 1.23 {dmA/dt)/{dV>U/dt) iiG 



(5) 



(c f. Equation [it and for Ve la is 22^G. 

Ivan Ommen et ah ( 1997 ) made a similar analysis for 
PSRs B1556-44 and B1727-47. For PSR B1556-44, 
the estimated electron density and the effective mean 
line-of-sight component of magnetic field in the vary- 
ing region are Ane = 290 cm~'^ and (-B||)var = 2 
directed away from the pulsar; for PSR B1727— 47, the 
electron density is Ang = 170 cm~^, the magnetic field 
strength is (-B||)var = 16 IJ-G directed toward the pul- 
sar. Because of the absence of an associated supernova 
remnant the authors proposed that the RM variation 
for PSR B1727— 47 may be due to a an irregularity 
within the Hii r egion of I Scorpius located at a dis- 
tance of 1.3 kpc. I Rankin et al. I (l988l) made a detailed 
analysis of the RM and DM variations for the Crab 
pulsar and found correlated variations between 1972 
and 1974 which imply a magnetic field in the vary- 
ing region of about 170 /xG. Several o ther authors (e.g. 
Weisberg et"al1l2004[ iHan et al.ll2006t ) also find appar- 



ent time variations in the interstellar RMs of several 
pulsars. 

Much stro nger fields are observed in vicinity of stars. 
For example, Johnston et al. ( 20051 ) observed RM vari- 



ations of many 1000 rad m ^ and DM variations of 
~ 20 cm"'^ pc for signals propagating through the disk 
of SS 2883, the companion star to PSR B1259-63, im- 
plying fields in the mG range. RM variations of a 
few 100 rad roT^ have been observed as t he path to 
PSR J1801-2304 passed close to the Sun (jOrd et al " 



20071 ) . With expected values of electron density in the 



solar corona, these variations again imply fields of order 
1 mC. 

The Parkes Pulsar Timing Array (PPTA) project is 
making regular timing observations of 20 millisecond 
pulsars (MSPs) with the primary goa l of detecting low- 
frequency (uHz) gravitational waves (|Manchesteiil2008 : 
iVerbiest et al.. .2010 ). Data from the first t hree years 



of thi s project (2004 - 2007) were used by lYou et al 



(|2007l ) to investigate the DM variations in the PPTA 
pulsars. Significant variations were found in many 
of the pulsars, with observed rates of change up to 



10 ^ pc cm 



yr ^. The PPTA sample is unusual 



if not unique in having a large number of observations 
over a long data span with high signal/noise ratio and 
accurately calibrated polarization data for a significant 
number of pulsars. Mean pulse total intensity and po- 
larization p rofiles have been published for the 20 PPTA 
pulsars by lYan et al.l (|2011l ) showing that the pulsed 



3 



emission extends over more than half of the pulse pe- 
riod for most of the sample and that the polarization 
variations across the profile are generally complex. 

In this paper we report on the long-term variations 
of observed mean PAs over a data span of 4.8 years 
for the 20 PPTA MSPs. It is shown that, after cor- 
recting for Faraday rotation in the Earth's ionosphere, 
there is little or no significant variation in the observed 
PAs and hence also in the interstellar RMs. Details of 
the observations and data processing are given in iJJl a 
comparison between two different ionospheric RM cor- 
rection methods is made in fj3l the observed variations 
are discussed in []4]and we provide a short summary in 



2 Observations and analysis 

The PPTA observations were made at intervals of 
2-3 weeks using the Parkes 64-m radio telescope 
with t he centre beam of the 2 0cm Multibeam re- 
ceiver |Stavelev-Smith et allEooj ). The Multibeam re- 
ceiver has a bandwith of about 300 MHz centered at 
about 1.4 GHz. The system equivalent fiux density on 
cold sky is approximately 30 Jy. The Multibeam re- 
ceiver has an orthogonal linearly polarized feed system. 
There is a calibration probe at 45° to the two signal 
probes through which a linearly polarized broadband 
and pulsed calibration signal can be injected. Three dif- 
ferent back-end systems were used, namely the Parkes 
digital filterbank systems PDFBl, PDFB2 and PDFB4. 
The PDFB systems are digital polyphase filterbanks us- 
ing field-programmable gate array (FPGA) processors. 
Signals from each of the two orthogonal polarizations 
are band-limited to 256 MHz and 8-bit sampled at the 
Nyquist rate. After transformation to the spectral do- 
main, typically with 1024 channels across the band, 
data are folded at the apparent pulse period to form 
mean pulse profiles in the four Stokes parameters. 

A ll data were record ed using the PSRFITS data for- 
mat ( Hotan et al .120041) with 1-minute sub-integrations. 
For PSRs J1857-h0943, J1939+2134 and J2124-3358, 
each observation was 32 minutes duration; for the other 
17 pulsars the observation times were 64 minutes. A 
short 2-minute observation of the pulsed calibration 
signal was made before each pulsar observation. Flux 
density scales were established using observations of 
Hydra A which was assumed to have a fiux density of 
43.1 Jy at 1400 MHz and a spectral index of -0.91. 
The Multibeam receiver has significant coupling be- 
tween the two signal feed probes. The cross-coupling 
was measured by analysis of series of observations of 
PSR J 0437— 4715 coverin g a wide range of parallactic 
angles (jvan Stratenll2004l) . 



All observation s were processed us ing the psrchive 
software package ( Hotan et al. 20041 ). Strong narrow- 
band radio frequency interference (RFI) was identified 
by searching for spectral channels that were signifi- 
cantly above the local mean; such channels were ex- 
cised from the data. Similarly, broad-band impulsive 
RFI was identified on zero-DM summed time-domain 
data and excised. The variations in instrumental gain 
and phase across the band and the effects of cross- 
couping in the feed were removed using the calibra- 
tion observations. Polarization parameters are in ac- 
co rdance with the astronom ical conventions described 
by van Straten et al. ( 2010l ): in particular, PAs are ab- 



solute and referred to celestial north, increasing toward 
east. Table[T]lists the basic parameters for the 20 PPTA 
MSPs with RM values derived from PDFB2 observa- 



tions following I Yan et al.l (|201lh . 

To measure the long-term changes in mean PA, we 
took the weighted average of the differences between 
the PAs across an individual observation profile and 
the PA s across the grand-average profile from lYan et al.l 
(|2011l ). This gives a relative PA shift for that observa- 
tion and its estimated uncertainty. This is then re- 
peated for all observations of a given pulsar to give the 
PA variations across the data span. In order to accu- 
rately align the observation profiles for a given pulsar, 
pulse times of arrival were obtained for all the obser- 
vations and fitted using th e TEMP02 pulsar timing 
package ijHobbs et al.ll2006l ) to give a timing model for 
the pulsar. This timing model was then used to align 
the profiles before computing the PA differences. 

Because the Earth's magnetic field is relatively 
strong (~ 0.5 G) and the electron density in the Earth's 
ionosphere relatively high (up to 10^ cm^"^), the Earth's 
ionosphere makes a signficant contribution to the total 
RM along the path to the pulsar. For the Southern 
Hemisphere, the Earth's magnetic field is directed out- 
wards and so the ionospheric contribution to the RM is 
negative, typically between —0.5 and —3 rad m~^, with 
a strong diurnal variation due to the varying ionspheric 
electron density. 



3 Ionospheric RM corrections 

We have used two programs to compute the ionospheric 
RM (RMi) in the direction of a given source at a given 
time. The first, parrot, was developed at the Domin- 
ion Radio Astrophysical Observatory (DRAG), Pentic- 
ton, Ganada. It uses three Ghapman layers to model 
the D, E and F regions and a semi-empirical model 
for variations in these based on the 10.7-cm solar radio 
flux measured by DRAG. The second, GEtRM_iono, 



4 



Table 1 Parameters for the 20 PPTA pulsars 
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was originally written by J. L 



Han for corr ection of 
2006[ ). In its 



Han et al 



measured pulsar RMs (e.g. 
current form, it uses the 2007 International Reference 
Ionosphere (IRI) modef^ with the URSI model for the 
critical frequency of the F2 ionospheric layer (/oF2) and 
input of the geomagnetic Ap index, the solar 10.7-cm 
radio flux and sunspot numbers. Both models use the 
International Geomagnetic Reference Fieltil to define 
the structure of the Earth's magnetic field. 

We have used these programs to compute RMi for 
each observation of each pulsar. Observed PA differ- 
ences have been corrected for the rotation occurring in 
the ionosphere to give the PAs above the ionosphere. In 
Figure [T]-|3] we show the observed PA differences in the 
upper panel for the 20 PPTA pulsars. The middle panel 
gives the PA differences after correction for ionospheric 
rotation using FARROT to compute RMi and the lower 
panel shows the PA differences after correction using 
getRM_iono. 

This figure shows that the PA variations are less af- 
ter removal of the ionospheric contribution. This is es- 
pecially noticable for pulsars in the northern sky, e.g., 



for 



^See http://iri.gsfc 
of the IRI. Source and data 



'■gov| 



files 



general description 
are available from 



|http://nssdcftp. gsfc.nasa.gov/models/ionospheric/iri/iri2007/] 
^ See [http://www.ngdc.noaa.gov/IAGA/vmod/igrf.htmll 



PSRs J1713-h0747 and J1939-H2134. For these pulsars 
viewed from Parkes the line of sight to the pulsar is 
roughly tangential to the local magnetic field (the dip 
angle at Parkes is about —60°) and so the RMi is large. 
There is a clear annual variation in the observed PAs 
resulting from the aliasing of the diurnal ionospheric 
variation with the annual sidereal motion. This annual 
variation is largely removed by the ionospheric correc- 
tion. Table [2] gives the rms scatter in the PAs and 
the corresponding reduced after a weighted fit of a 
straight line to the corrected PAs for both the farrot 
and IRI corrections for the PPTA pulsars. Without ex- 
ception the PAs after the IRI corrections have a smaller 
rms scatter and reduced than those after the FAR- 
ROT correction. 



4 Results and Discussion 

After the ionospheric correction, the PA variations for 
essentially all of the PPTA pulsars are dominated by 
a short-term scatter. As shown by the large values of 
reduced (Table [2]) , this scatter is much larger than 
the nominal uncertainty of the PAs in most cases. For 
the stronger pulsars, the PA uncertainty estimated from 
the off-pulse rms variations in the Stokes profiles is very 
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Fig. 1 PA differences as a function of time for six of the PPTA pulsars. In eacii sub-figure, the upper panel shows PA 
differences without ionospheric correction, the middle panel gives PA differences after ionospheric corrections using the 
Penticton program PARROT; the lower part shows differences after ionospheric corrections using the IRI ionospheric model. 
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Fig. 2 PA difFerences as a function of time for six of the PPTA pulsars. See Figure [T] for details. 



7 





• • • 


• • • • 

• • 
• 


3 o 
< 

Ph 


• + 


•« • 


1 1 1 1 1 1 1 1 1 1 1 

.+ 

•• • Jb* 


1 1 1 1 1 1 1 1 1 1 
i 

• • « 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 

• •*•%••••* 


S o 
< 

Ph 


1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 

+ 

J1744-1134 


1 1 1 1 1 1 1 1 1 1 
• • 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


a, 
< 


1 1 1 1 1 1 1 1 1 1 

Jia34-345a 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


-750 -500 


-250 

MJD-54550 


250 500 750 




-750 -500 


-250 250 500 750 
MJD-54550 




• * 

W 


.+ . A ++.•!, + _ 

+ - 


< 

Ph 






1 1 1 1 1 1 1 1 1 1 1 

- 

+ + 


1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 
+ - 


S o 
< 

Ph 


1 1 1 1 1 1 1 1 1 1 

±^ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

!• : • • 


1 1 1 1 1 1 1 1 1 1 1 

- + ^^v* 

Jia57+0943 


1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

+ - 


a, 
< 


1 1 1 1 1 1 1 1 1 1 

J1909-3744 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

• • • • * 

r 


-750 -500 


-250 

MJD-54550 


250 500 750 




-750 -500 


-250 250 500 750 
MJD-54550 






• "'l 
• . • ••• • • 


< 

Ph 




T 


1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


in 

3 o 
< 

Ph 


1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 

J1939+3134 
-750 -500 


1 1 1 1 1 1 1 1 

-250 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

* s * 

250 500 750 


a, 
< 


1 1 1 1 1 1 1 1 

J2124-335B 
-750 -500 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



Fig. 3 PA difFerences as a function of time for six of the PPTA pulsars. See Figure [T] for details. 
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Fig. 4 PA differences as a function of time for two of tlie PPTA pulsars. See Figure [T] for details. 



Table 2 Rms variations and reduced for ionosphere-corrected PAs 



PSR FARROT IRI 





CPA 

n 


Reduced 

x' 


OTA 

n 


Reduced 

x' 


J0437-4715 


1.14 


833.3 


0.96 


488.6 


J0613-0200 


2.31 


5.9 


2.15 


4.7 


J0711-6830 


1.82 


3.5 


1.70 


3.3 


J1022+1001 


1.22 


192.8 


0.93 


96.5 


J1024-0719 


1.41 


47.3 


1.31 


43.1 


J1045-4509 


1.74 


5.6 


1.63 


4.6 


J1600-3053 


1.33 


10.6 


1.08 


6.5 


J1603-7202 


1.35 


4.0 


1.29 


2.9 


J1643-1224 


1.43 


5.7 


1.09 


3.5 


J1713+0747 


1.40 


76.9 


1.37 


70.7 


J1730-2304 


1.83 


25.2 


1.71 


17.0 


J1732-5049 


1.93 


1.8 


1.92 


1.5 


J1744-1134 


1.28 


779.9 


0.98 


411.9 


J1824-2452 


1.08 


16.5 


0.95 


11.1 


J1857+0943 


3.02 


5.3 


2.71 


4.0 


J1909-3744 


1.15 


138.1 


1.11 


123.5 


J1939+2134 


1.22 


64.3 


1.42 


50.6 


J2124-3358 


1.44 


3.4 


1.31 


3.0 


J2129-5721 


1.72 


5.4 


1.71 


5.2 


J2145-0750 


2.49 


31.8 


2.31 


20.2 



9 



small, often ^ 1°. Residual calibration errors can ex- 
ceed the rms noise in the stronger pulsars, leading to 
PA offsets and hence scatter in the PAs. Low-level un- 
excised RFI can affect both the accuracy of the cali- 
bration and the Stokes data themselves also leading to 
PA scatter. Finally, although the ionospheric correc- 
tion has reduced the short-term scatter in the PAs, it 
cannot be completely accurate as it relies on models 
for the ionosphere and its variations whereas the actual 
ionosphere is dynamic and changes on short timescales 
in often unpredictable ways. It is difficult to separate 
these contributions to the scatter in the corrected PAs, 
but it is probable that this last effect is dominant. Near 
real-time monitoring of the total electron content, for 
example, usin g the dual-frequenc y Global Positioning 
System (e.g., Garner et al. 20081) . should help to im- 
prove the accuracy of the RMi measurements and hence 
reduce the overall short-term scatter in the PAs. 

Figures [T] - 2] show that there is little if any signif- 
icant long-term variation in the corrected PAs. This 
implies that the mean PA of the emitted pulses and 
hence the orientation of the spin and magnetic axes of 
the pulsar have varied hy ^ 1° over the nearly 5-year 
data span. This is not surprising since six of the 20 
pulsars are isolated and none of the 14 binary pulsars 
is in a clo se orbit where geodetic precession might b e 
expected ( Breton et al. 2008 : Manchester et al. 2010l ). 
In addition, most of the binary companions are evi- 
dently of relatively low mass (J; 0.5 Mq). Because of 
the great age of these MSPs, any free precession result- 
ing from interactions or deformations in the spin-up 
phase would be likely to have dissipated by now (cf. 
Glampedakis et al. 20091) . 

Figures [5] - [8] show the IRI-corrected PAs for the 
PPTA pulsars. The straight line on each plot results 
from a weighted least-squares fit of a straight line to 
the data. To reduce the short-term scatter, weighted 
averages of the corrected PAs over 100-day blocks are 
shown in the middle panel, again with the fitted straight 
line. If we assume that there has been no variation in 
the intrinsic PA over this interval, we can convert the 
observed PA variations to equivalent RM variations us- 
ing Equation [T] as shown in the lower panel. Table |3] 
gives the slope of the fitted line and its estimated un- 
certainty and the corresponding RM gradient for the 
20 pulsars. Only a few of these gradients are significant 
at more than the 3cr level; mostly the observations are 
consistent with no long-term change in PA or interstel- 
lar RM within the uncertainties. 

As discussed in ^Jl] we can use the ratio of the RM 
gradient to the DM gradient to estimate the line-of- 
sight component of the magnetic field in a putative 
time- varying region. D M variations for the PPTA pul- 
sars were measured bv lYou et al.l (|2007| ). There is only 



one year of DM data overlap with the PA variation 
data but the DM gradients given in Table [Sj obtained 
by fitting to this overlap year, are representative of the 
longer-term gradients. Because of the high-precision 
timing possible for these MSPs, the absolute changes 
in DM, while significant in some cases, are extremely 
small. This, combined with the fact that we effectively 
have just upper limits on dDM/dt for most of our sam- 
ple, means that, for these pulsars, we can just place 
upper limits of several 100 /zG on the interstellar field 
strength in the time-varying region. Values and uncer- 
tainties of (i3||)var given in Table[3]have been computed 
using Equation [S] with the uncertainties in dDM/dt and 
dKM/dt assumed to be independent. 

Apparently significant long-term PA and implied 
RM variations (> 3cr) are observed for five of the 
pulsars: PSRs J0613-0200, J1022+1001, J1603-7202, 
J1909-3744 and J2145-0750. For PSRs J1022+1001 
and J2145— 0750, dDM/dt is not well determined and so 
the mean field is very uncertain. For PSRs J0613— 0200, 
J1603-7202 and J1909-3744, there are significant val- 
ues for both dKM/dt and dDM/dt, so we can compute 
values for (i?||)var as listed in Table [3] For most of the 
pulsars, the (i?||)var values or limits are much larger 



than typical interstellar fields of a few /zG (|Han et al 



|2006[) and even the field expected in ma gnetic filaments 
associated with superno va remnants (|Hamilton et al 
19851 iRankin et al.lll988l) . 



For PSRs J1024-0719, J1045-4509 and J1824-4509 
the upper limits on (i?||}var are interestingly low, of 
order 10 /iG. This is comparable to the estimated 
interstellar total field strength, including both large- 
scale regular fields ari d small-scale irr egular fields 
( Beuermann et al.l [l985t iHan et al.l |2004| ). These re- 
sults show that variations in the interstellar RM for 
these pulsars are not dominated by any strong filament 
along the line of sight. 

As well as the significant long-term variations in RM 
discussed above, Figures [S] - [5] show that there are ap- 
parently significant fluctuations in the smoothed RM 
variations on a timescale of a few hundred days for some 
of the MSPs. For PSR J1939-h2134, the PA appears to 
increase at a rate of about 3 rad m~^ yr^^ for the first 
150 days (MJD 53900 - 54050). Over this same interval, 
the DM is decrea sing at about —8 x 10"'* cm~'^pc yr~^ 
( You et all |2007| ). If both these variations are occur- 
ring in a single region, they imply a line-of-sight mag- 
netic field of about 4.5 mG in the region. Although 
such fields are observed close to stars, such an inter- 
section is highly improbable (except for the Sun - dis- 
cussed below). Most probably, the observed variations 
are not common to a single region, but represent statis- 
tical fluctuations in the accumulated RM and DM over 
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Fig. 5 PA variations for six of tlie PPTA pulsars. The upper panel of each sub-figure has the IRI-corrected PA variations 
showing the fitted linear trend. The middle panel shows the same data averaged over 100-day blocks with the fitted line as 
in the upper panel. The lower panel shows the smoothed variations converted to RM changes using Equation ([l]). 
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Fig. 6 PA variations for six of the PPTA pulsars. See Figure [5] for details. 
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Fig. 7 PA variations for six of the PPTA pulsars. See Figure [5] for details. 
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Fig. 8 PA variations for two of the PPTA pulsars. See Figure [5] for details. 



Table 3 PA gradients, implied RM gradients, DM gradients and apparent line-of- 
sight magnetic field in the time- varying region for the PPTA pulsars. 



PSR 


dPA/dt 


dRM/dt 
(rad m~^ yr "'^) 


dDM/dt 
(cm^^ pc yr^^) 


(-B||)var 

(mG) 


J0437-4715 
J0613-0200 
J0711-6830 
J1022+1001 
J1024-0719 


0.08 ± 0.03 
0.64 ±0.13 

-0.08 ±0.13 
0.35 ±0.06 

-0.08 ± 0.09 


0.03 ±0.01 
0.23 ±0.05 

-0.03 ±0.05 
0.13 ±0.02 

-0.03 ±0.03 


(3 ± 3) X lO-"^ 
(4.3 ± 1.6) X 10"-* 
(-5.3 ±7.8) X lO""* 
(-2.0 ±5.5) X 10-"* 
(4.8 ± 1.9) X 10"^ 


1200 ± 1300 
660 ± 280 
70 ± 160 
-800 ± 2200 
-8 ±8 


J1045-4509 
J1600-3053 
J1603-7202 
J1643-1224 
J1713+0747 


-0.12 ±0.09 
0.05 ±0.07 
0.33 ± 0.08 
-0.11 ±0.08 
-0.08 ±0.08 


-0.04 ±0.03 
0.02 ±0.02 
0.12 ±0.03 
-0.04 ±0.03 
-0.03 ±0.03 


(-1.03 ±0.01) X 10"^ 
(-9± 10) X 10"" 
(3.4 ±0.8) X 10-3 
(-3.2 ±3.7) X 10"* 
(-6 ±24) X lO"'^ 


5±4 
-30 ± 40 
40 ± 10 
150 ±210 
620 ± 2500 


J1730-2304 
J1732-5049 
J1744-1134 
J1824-2452 
J1857+0943 


0.17 ±0.12 
0.47 ±0.30 
0.04 ± 0.06 
-0.02 ±0.10 
-0.35 ±0.21 


0.06 ±0.05 
0.17 ±0.11 
0.01 ±0.02 
-0.01 ±0.04 
-0.13 ±0.08 


(-2.3 ±7.1) X lO""* 

(7± 11) X 10"* 
(-5.0 ±8.7) X lO""* 
(6.2 ±0.7) X 10-=* 
(-1.4 ±3.3) X 10-3 


-320 ± 1000 
300 ± 510 
-20 ± 70 

-2 ±8 
110 ± 280 


J1909-3744 
J1939+2134 
J2124-3358 
J2129-5721 
J2145-0750 


0.37 ±0.07 
-0.25 ± 0.11 
-0.04 ± 0.07 

0.25 ±0.10 
-0.78 ±0.14 


0.13 ±0.02 
-0.09 ±0.04 
-0.02 ±0.02 

0.09 ± 0.04 
-0.28 ±0.05 


(-4.6 ±0.9) X 10-* 
(3.3 ± 1.7) X 10-* 

(-4.5 ±8.1) X 10-* 

(5 ±4) X 10-* 
(3.5 ±4.6) X 10-* 


-350 ± 90 
-340 ± 230 
50 ± 110 
220 ± 200 
-990 ± 1300 
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the whole path resuhing from random and uncorrelated 
spatial and temporal variations in the interstellar elec- 
tron density and magnetic field. 

It is important to note that the observed RM varia- 
tions would be expected to be less than the DM vari- 
ations since the RM is a vector sum (Equation ^) 
whereas the DM is a scalar sum (Equation Vari- 
ations in regions of oppositely directed field will cancel 
in RM but not in DM. A more detailed examination 
of the origin of these variations would require realis- 
tic modelling of the electron density and magnetic field 
fluctuations along representative Galactic paths. This 
is left for future work. 

It is also worth noting that the aliased diurnal vari- 
ations in ionospheric RM are approximately covariant 
with annual variations in RM and DM resulting from 
the changing path through the solar corona as the Earth 
moves around the Sun. Both have maxima when the 
pulsar transits near local noon. At least for paths that 
do not pass too close to the Sun, the DM variations 
are roughly sinusoidal; Tempo2 has a correction term 
for d e lays occurring in the solar corona ([Edwards et al. 
2006t I You et al.llioOTal ). Predicted RM variations are 
more com plex because of structure in the solar mag- 
netic field |Ord et ahlboOTl ). Some of the PPTA pulsars 
(notably PSRs J1022-hl001 and J1730-2304) have low 
ecliptic latitudes, but the few observations when these 
pulsars were close to the Sun have been excluded from 
the present analysis. The fact that little or no annual 
PA variation is evident in Figures [S]- [5] shows that any 
variations due to the changing path through the solar 
corona are below the level of detectabilty for these ob- 
servations of the PPTA pulsars. 

The main goals of the PPTA and other similar PTA 
projects depend on precise measurement of pulse times 
of arrival (ToAs). In principle, the ionospheric RM 
should be taken into account when summing the data 
in frequency. Failure to do so would result in some 
depolarization of the signal. Depending on the ToA 
algorithm used and the details of the profile polariza- 
tion this could affect the ToAs. Fortunately, at the 
frequencies of 1.4 GHz and above commonly used for 
PTA observations, the effect is very small, with the de- 
polarization being 1 per cent across the band. Some 
To A algorithms, for e xample, matrix template match- 
ing ( van Straten 20061) . solve for the PA rotation across 
the band and so the effect is decoupled from the ToA 
measurement. This is potentially important for pulsar 
timing at low radio frequencies. 



5 Conclusions 



for investigation of long-term variations in the proper- 
ties of the observed pulses. In this paper we have anal- 
ysed the time variations of mean absolute pulse position 
angle over 4.8 years. The largest systematic effect found 
was that due to variations in the Earth's ionosphere. 
Two ionospheric models, farrot and the International 
Reference Ionosphere (IRI) model, were used to predict 
the variations in the ionospheric component, RMi, and 
to compensate the observed PA changes for ionospheric 
variations. It was found that the IRI model gave better 
results than farrot. 

After compensation, the rms fluctuation in the PAs 
was about 1° with little or no signiflcant long-term vari- 
ations over the data span for most pulsars. This implies 
not only that the polarization of the emitted radiation 
was very stable, but also that the net long-term vari- 
ations in interstellar RM along the path to the pulsar 
were also very small, typically less than 0.1 rad m^^ 
yr~^ in absolute value. The observed variations were 
dominated by short-term fluctuations that most likely 
originate from residual calibration errors or errors in 
the RMi due to unmodelled short-timescale variations 
in the ionospheric total electron content. 

Apparently significant PA/RM variations on timescales| 
of a few 100 days or longer were observed in some of 
the pulsa rs. Comparison w ith the DM variations pub- 
lished bv lYou et all |2007[) showed that these RM vari- 
ations imply high values of magnetic field if they are 
interpreted as correlated variations occurring in a small 
part of the path. More likely, they represent statisti- 
cal fiuctuations in the whole-of-path RM resulting from 
random spatial and temporal variations in both the in- 
terstellar electron density and magnetic field. 
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The regular observations of 20 MSPs undertaken as 
part of the PPTA project provide an excellent data base 
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